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Cooper pairing of spatially separated electrons and holes in graphene bilayer is studied beyond the 
mean-field approximation. Suppression of the screening at large distances, caused by appearance 
of the gap, is considered self-consistently. A mutual positive feedback between appearance of the 
gap and enlargement of the interaction leads to a sharp transition to correlated state with greatly 
increased gap above some critical value of the coupling strength. At coupling strength below the 
critical, this correlation effect increases the gap approximately by a factor of two. The maximal 
coupling strength achievable in experiments is close to the critical value. This indicated importance 
of correlation effects in closely-spaced graphene bilayers at weak substrate dielectric screening. An- 
other effect beyond mean-field approximation considered is an influence of vertex corrections on the 
pairing, which is shown to be very weak. 

PACS numbers: 73.22.Pr, 71.35.-y, 74.20.-z, 73.21.-b 



I. INTRODUCTION 



Cooper pairing of spatially separated electrons and 
holes due to their Coulomb attraction was proposed ini- 
tially as apossible origin of superfluidity [IJ, Josephson 
effects [ll-y] and anomalous electromagnetic phenomena 
0, [1] in coupled semiconductor quantum wells. Strong 
Coulomb electron-hole attraction was supposed to main- 
tain high critical temperature of the pairing, while spa- 
tial separation of paired electrons and holes could prevent 
them from interlayer tunneling leading to their recombi- 
nation and to a condensate phase fixation. Experimen- 
tal evidences of superfluid transition in two-layer semi- 
conductor structures in strong magnetic field (pairing of 
composite fcrmions) 0], and without magnetic field 7- 
Il2l | were found several decades later. 

Experimental fabrication of graphene, an atomically- 
thin two-dimensional form of carbon [isl . [T^ . opens a 
possibility to realize electron-hole pair condensation in 
a spatially separated graphene bilayer. Structures con- 
sisting of two independently gated graphene layers with 
common contacts and very small (0.6 nm) separation, 
have been fabricated and studied experimentally fisl [l6l | . 
A system of two independently contacted graphene lay- 
ers separated by 5 nm-thick Si02 barrier have been also 
made [l^l- Coulomb drag have been already studied 
in such system. The most promising for realization of 
the pairing could be heterostructures consisting of two 
graphene layers separated by atomic-thin boron nitride 
layer. Considerable progress in fabrication of such struc- 
tures with ultra-high mobility encapsulated graphene 
samples have been achieved recently [1&~2^] . Measure- 
ments of Coulomb drag in these structures have been 
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reported (see [2j| and references therein). 

Theoretical studies of electron-hole Cooper pairing in 
spatially separated graphene bilayer were presented in 
|24l - l33| . It was shown that both weak and strong coupling 
regimes of the pairing are achievable experimentally . 
Estimates of the critical temperature Tc at strong cou- 
pling [25 'SC?! are very different depending on model and 
approximations used: ranging from room temperature, 
with the unscreened Coulomb attraction taken as a pair- 
ing potential [IE Hsj , to unobservably small values given 
by Bardeen-Cooper-Shrieffer (BCS) theory with screened 
interaction p7{ . 

It was argued in [H, [H, [s^ that the pairing in 
graphene bilayer at strong coupling is multi-band, i.e. 
involving both conduction and valence bands of elec- 
trons and holes. The reason is that electrons and holes 
in graphene are described by effective two-dimensional 
Dirac-type equation for massless particles [l^, ^-nd 
thus a gap between conduction and valence band is ab- 
sent. It was shown [2^, [2^ [s^ that in this case can 
be much larger than given by one-band BCS-like model 
(27j . The study of the problem with taking into account 
frequency-dependence effects was presented in [30| . 

Mean-field approximation, used in most of these theo- 
retical works [24 30, 33], is known to be well- applicable 
at both weak- and strong-coupling sides of BCS-BEC 
crossover, taking place in conventional pairing systems 
[H,!!!!. However, little is known about its applicability 
in the regime of multi-band Cooper pairing, which takes 
place in graphene bilayer at strong coupling instead of 
the BEC (Bose- Einstein condensate) side of the crossover 
[29I l36j (the BEC regime restores iri graphene bilayer if 
a gap is opened in its spectrum [37L [38i|). 

One of correlation effects arising beyond the mean-field 
approximation is suppression of the screening (and con- 
sequent enlargement of electron-hole interaction) due to 
appearance of the gap and order parameter in the sys- 
tem. As supposed in 3l|, mutual positive feedback be- 
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tween onset of the gap and suppression of the screening 
can result in appearance of two solutions of the gap equa- 
tion: small gap at strong screening and large gap at weak 
screening. 

Change of polarizability of electron system caused by 
appearance of a gap in its spectrum have been exten- 
sively studied in a context of collective modes in super- 
conductors (see, e.g., [3^14^]) and phonon self-energies, 
acquiring sharp features at frequencies twice larger than 
the gap (see [SI 113 references therein). 

In systems with electron-hole pairing (proposed ini- 
tially as an origin of excitonic insulator [43| ) suppression 
of a system polarizability due to the pairing can change 
significantly the pairing interaction itself. The question 
of self-consistent treatment of pairing and screening in 
electron-hole liquids in semiconductors have been ad- 
dressed in several works [43 - l48| . 

For electron-hole graphene bilayer, the self-consistent 
suppression of the screening, suggested in [3l|, was con- 
sidered in (32I with using the one-band approximation, 
however, without details of calculations presented. The 
result was that this effect is negligible for a considered 
range of parameters and the gap equation has always 
only one solution, corresponding to a very small gap. 
In our article, we study this effect within the model of 
multi-band pairing at strong coupling. We show that (in 
contrast to [S^l) this effect can drastically change charac- 
teristics of the pairing at strong enough coupling, leading 
to formation of very large gap. 

The self-consistent suppression of the screening is de- 
scribed by a series of Feynman diagrams, involving 
the screening of electron-hole interaction by virtual Bo- 
golyubov excitations. Considerable part of remaining 
diagrams beyond mean-field Gor'kov equations can be 
absorbed into renormalization of a Coulomb interaction 
vertex. In the theory of superconductivity this renor- 
malization is negligible at weak coupling due to Midgal 
theorem, but can be appreciable at strong coupling, gen- 
erally increasing a critical temperature [49|. Corrections 
to Coulomb interaction vertex in graphene were consid- 
ered in |50H53| . It was shown that the vertex corrections 
are large in undoped graphene 5l|, but rather small at 
finite doping [52I]. We calculate numerically the simplest 
vertex correction and show that it is small, though acts 
to increase a coupling strength. 

The article is organized as follows. In Sec. 2 we intro- 
duce briefly the multi-band model of the pairing, in Sec. 3 
we consider self-consistent suppression of the screening 
and its effect on the pairing at zero temperature. Sec. 4 
presents calculations of the vertex corrections and Sec. 5 
concludes the article. 



[28l I29I |33| . Here we present only the formulas needed 
for further calculations. 

To describe the pairing, we introduce the Matsub- 



ara Green functions (p, t) — — (rcp-|^(r)cp^2*"(0)), 

where Cp^ = Op^, Cp^ = ap^-j, Op-y and ap^ are destruc- 
tion operators of electrons in electron- and hole-doped 
layers respectively with the momentum p from conduc- 
tion (7 = -1-1) or valence (7 = —1) band. The bare Green 
functions for electrons and holes are: G^|,"72 (Pj *£«) — 



where = -^^^^ = $p 



P7 - -S,P7 = S,P7 - - M 

are energies of electrons and holes measured from the 
chemical potentials fi and —fi in electron- and hole-doped 
graphene layers respectively; vp « lO^m/s is the Fermi 
velocity in graphene. 

In the simplest case of s-wave band-diagonal pairing, 
all the Green functions G^^l^^ are nonzero only at 71 = 72 ■ 
Solution of the Gor'kov equations in the Cooper channel 
in this case is: 
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and G^^y^' (p, ie 



G^jj'{p,ien)- Conventional notations for coherence fac- 
tors and Bogolubov energies are used: 
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Note that in the multi-band regime we have two kinds of 
the gap functions Ap± and Bogolyubov excitation ener- 
gies Ep± corresponding to conduction and valence bands. 
The gap functions are determined by equations: 
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The factor F^J, ^ |(p7|p'7')l' = (1 + u'pp' /pp')/2 
arises in result of summation over spinor components of 
the effective electron wave function in graphene [1^ . We 
use here the static approximation, in which the gap Apj 
and the statically screened interaction V{q) are real and 
independent of frequency. Substituting ([I])-© into 
performing summation over e„ = 7rT(2n -I- 1) and inte- 
gration over an angle of p', we get at T — >• 0: 



II. MULTI-BAND DESCRIPTION OF THE 
PAIRING 



Apy — 



7 



00 



A 



p'7' 



2 En 



(4) 



Multi-band pairing of electrons and holes in graphene 
bilayer at strong coupling is described in detail elsewhere 



where U^^Mp,/) = [Voip,p') + -fj'Vi{p,p')]/2 is a half- 
sum of half-difference of s- and p-wave harmonics of the 



3 



potential: 



Vi{p,p') ^ 1 Y' cos(Z'^)F (^^Jp^ + - 2pp' cos ipj . (5) 



As shown in |29|] , the first harmonic Vi can be neglected 
at strong coupling. In this case the gap functions in con- 
duction and valence bands are equal. We will find an ap- 
proximate solution of (U in a manner similar to that used 
in [33|: we assume that Apj = A x f{p), where f{p) is 
some trial function, satisfying the conditions /(pf) = 1, 
f{p)(x 1/p at p — > oo. Fixing in ^ p = pp and neglect- 
ing Vi, we get the following algebraic equation for the 
gap A at r = 0: 



n/N 





1 

- 

— 


I \ \ \ 


- 








- 










1 , 1 





1 = 



p'dp' 
Stt 



Vo{pF,p')fip') 



1 



1 



(6) 



III. SELF-CONSISTENT SUPPRESSION OF 
THE SCREENING 

The statically screened potential of electron-hole in- 
teraction, entering ([S]), in random phase approximation 
(well-applicable in graphene bilayer due to large degree 
of electron state degeneracy [l^, IH, 'Hi] ) is [23. l33j| : 
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where Vq = jeq is the bare Coulomb attraction, 
screened by surrounding medium with dielectric per- 
mittivity e, no((j') is a static polarizability of a single 
graphene layer, D is the interlayer distance. 

The most favorable conditions for the pairing are 
achieved at small interlayer separation, when p^D <C 1. 
In this case a behavior of the system is determined only 
by the dimensionless parameter = e^/ewp ~ 2.19/e, 
which specifies the coupling strength, while ([7]) reduces 
to 



^9) = 
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The static polarization operator of doped 
graphene, calculated without taking into ac- 
count electron- hole pairing, is [ssl ISa ] no(g) = 
5AA{-l-He(g-2pF)(<z/4pF)G<(2pF/g)}. Here g = 4 
is the degeneracy factor, M = /i/27rw| is the density of 
states (per spin projection and valley) at the Fermi level, 
G<{x) — xy/1 — — arccosx. The long- wavelength 
asymptotics 110(9) ~ provides metallic-like 

screening at long distances. 

Formation of a condensate of interlayer electron-hole 
pairs leads to appearance of a direct response of charge 
density in one layer on electric field in the other layer, 
described by anomalous polarizability Ha. Normal, in- 
tralayer polarizabilities n,i of each graphene layer also 



FIG. 1: (Color online) The effective static polarizability 11 
(in units of A/") in graphene bilayer with Cooper pairing as a 
function of momentum q at different values of A//^, indicated 
above corresponding curves (solid lines). Dotted line: the 
intrinsic static polarizability Ho at A = 0. 



change with respect to intrinsic polarizabilities IIq due to 
appearance of the gap in the energy spectrum. The ex- 
pression for y(g) at p^D <C 1 takes the same form as ([5]), 
but with the replacement no((?) — )■ n((7) = X\^{q)+X\^((j). 

In the random phase approximation the normal Iln 
and anomalous Ha polarizabilities are calculated as loops 
consisting of two normal or anomalous Green functions 
respectively: 
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where p' = p + q. Substituting ([T]) into ©-((101), per- 
forming frequency summations and taking a limit T — 0, 
we obtain the following expressions for the static polar- 
izabilities n„^a(g) = IIn^a(9, iw„ 0): 
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The sum 11 = !!„ -I- Ha, playing the role of effective 
polarizability, is plotted in Fig. [1] as calculated numeri- 
cally according to PT|) -([T ^ at different values of A. It is 
seen that polarizability of the system with Cooper pair- 
ing is suppressed at small momenta due to appearance 
of the gap. A magnitude and momentum region of this 
suppression grow at increasing A. The long- wavelength 
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FIG. 2: (Color online) Solid lines: the right-hand side I of the 
gap equation ([G]) /(A) = 1 as a function of A (in logarithmic 
scale) at different values of rs, indicated in the right. 

asymptotics Ti-{q) ~ —gM jVl'p^K^ indicates that the 
screening at long distances is absent when A 7^ and, 
moreover, transition from n((7) to no((7) at A is not 
uniformly continuous. 

One of the consequences of this discontinuity is 
that the usual BCS-like recipe, involving replacement 
Vo(p,p') — > K)(pf,Pf) in the gap equation (g]) becomes 
inapplicable, since the Fermi-surface value Vo(pf,Pf) di- 
verges due to absence of a long-range screening. However, 
the method ^ of approximated solving of the gap equa- 
tion is applicable even in this case, since the singularity of 
Vb(pF,p') at y' — pf is logarithmic and thus integrable. 

In numerical calculations, we use the trial func- 
tion f{p) — 1/{\p/pp — 1\ + 1) and the following ap- 
proximation for the effective polarizability: n(g) « 
-gAfq^ / [{12 A'^/viy^^ + q'^/^]^/^, which is close to the 
numerically calculated Tl{q) at any A and retains its ma- 
jor features — the correct asymptotic at smallest q and 
tenency to —gAf at larger q. 

The right-hand side / of Eq. ^ is plotted as func- 
tion of A at various in Fig. [2j The points of inter- 
section /(A) = 1 give solutions A of the gap equation. 
In the mean-field approximation, the dependence /(A) 
should be monotonously decreasing. Correlation effects 
make /(A) non-monotonous at rg ^ 0.2. This can po- 
tentially result in appearance of three solutions of the 
gap equation. Among these solutions, only the largest 
A, corresponding to the lowest ground-state energy, will 
be established in the system. 

The maximal gap A is plotted in Fig. |3] as a function 
of Tg . At small enough , the gap value is approximately 
twice larger than without taking into account correla- 
tion effects. When rg exceeds some critical value (about 
2.35 in our case), three solutions of ([6]) appear, and the 
maximal gap becomes very large and comparable to the 



FIG. 3: (Color onhne) Sohd hue: the largest value of the gap 
A (in logarithmic scale), found from Eq. ((6| at different values 
of Ts. At Ta ~ 2.35 the gap jumps (dashed line) to much larger 
value due to appearance of three solutions of Eq. Dotted 
line: the gap, calculated in the mean-field approximation. 

chemical potential. The critical value of depends on 
details of the model used (in particular, on a form of the 
trial function /(p)). In our case it is rather close to the 
maximal ~ 2.19 achievable experimentally at £ = 1, 
when graphene bilayer is suspended in vacuum. However, 
correlation effects can play important role in the vicinity 
of transition to the strongly-correlated state even for Vg 
smaller than the critical value. 



IV. VERTEX CORRECTIONS 

The bare vertex of electron-electron Coulomb interac- 
tion in graphene (Fig. [Ja)), entering the Gor'kov equa- 
tions, is ry^^(p',p) — (p'7'|p7), where the graphene- 
specific angular factor is 29] 

[ zsm^*^, 7 = -7 , 



(a) (b) 




Pi p 



FIG. 4: The bare vertex (a) of electron-electron Coulomb 
interaction and the simplest correction (b) to it. 
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the effect of vertex corrections on the pairing, we compare 
the second-order correction to the couphng constant on 
the Fermi surface — J\fU^l{pF , Pf) with its first- 

order value A^|,°l = 7VC/^°!(pf,Pf)- 

In Fig. [S] the conduction-band component of the cor- 
rection a'i^^ is plotted as function of rg. Its other com- 
ponents (aL^I and A^^^L^) are very close to it and thus 
are not shown. For comparison, the first-order coupling 
constants in the intra-band (A^^^) and inter-band (A^^L-^) 
channels are also plotted (these quantities were studied 
in detail in [29]). It is seen that the second-order vertex 
corrections amount to about 5% of the first-order cou- 
pling constants and thus can be neglected. 



FIG. 5: (Color online) The dimensionless second-order correc- 
tion A^^:^ (dotted line) to the coupling constant as a function 
of rs in comparison with the first-order intraband X'^j (solid 
line) and interband xS^''-.y (dashed line) coupling constants. 

ip and </j' are the azimuthal angles of p and p'. 

The simplest correction F*^^^ to F*^"^ is shown diagram- 
matically in Fig. IDJb). We take the bare Green func- 
tions (p2, je2) and G^'^^j (pi, iei) as internal elec- 
tron lines of this vertex, and the Coulomb interaction as 
the internal wavy line. Moreover, we consider only the 
static vertex, i.e. the vertex function at zero frequen- 
cies: F^^^(p',p) = F^^^(p',ze' ^ 0,p,ie ^ 0). After 
frequency summation in the internal loop, we get in a 
zero-temperature limit T — > 0: 

^ylyip'^P) = Y1 I 7^^(pVIP272>(P272|Pi7i) 

x(pi7i|P7>niPi-p|)^%^^^|^, (13) 

SP171 SP272 

where p2 = Pi + p' — P and Q{x) is a unit step function. 

The integral in diverges at small q if we take the 
bare Coulomb interaction V{q) — /eq. Thus we will 
take the statically screened interaction ([5]) as V{q). 

Substitution of the vertex correction (fT3|) into Gor'kov 
equations results in the correction to the potential 

u^%,p') = J ^(P7|pV)r5^(p',P), (14) 


added to u!"^2{p,p') in the gap equation (U). To estimate 



V. CONCLUSIONS 

We have considered the correlation effects in Cooper 
pairing of spatially separated electrons and holes in 
graphene bilayer at strong coupling. The first effect con- 
sidered is the self-consistent suppression of the screening 
of electron-hole interaction due to appearance of the gap 
and order parameter. Its most remarkable consequence is 
an absence of screening at long distances at any nonzero 
gap, that makes a usual BCS method inapplicable due to 
divergence of the interaction at the Fermi surface. 

Thus we performed numerical calculations with full 
momentum integration to solve the gap equation. We 
found that, at small enough coupling strengths, correla- 
tions increase the gap by a factor of two. However, at 
coupling strength above some threshold the gap sharply 
increases by several orders of magnitude, which indicates 
transition of the system to a strongly-correlated state. 
The critical value of the coupling strength depends on 
details of the theoretical model. In our case it turned 
out to be only slightly larger than the maximal value 
achievable in experiments. 

We have also considered the role of vertex corrections 
at strong coupling. The simplest correction to the in- 
teraction vertex was calculated numerically in the static 
approximation and was shown to enhance the pairing, 
but only on about 5%. 
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